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We present an exact solution for the distribution of sample averaged monomer to monomer dis¬ 
tance of ring polymers. For non-interacting and weakly-interacting models these distributions cor¬ 
respond to the distribution of the area under the reflected Bessel bridge and the Bessel excursion 
respectively, and are shown to be identical in dimension d> 2. A symmetry of the problem reveals 
that dimension d and A — d are equivalent, thus the celebrated Airy distribution describing the areal 
distribution of the d = 1 Brownian excursion describes also a polymer in three dimensions. For a 
self-avoiding polymer in dimension d we find numerically that the fluctuations of the scaled averaged 
distance are nearly identical in dimension d = 2, 3 and are well described to a first approximation 
by the non-interacting excursion model in dimension 5. 
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The statistical mechanics of polymers has been well 
studied for many years due both to the numerous practi¬ 
cal applications of polymers as well as their many inter¬ 
esting properties. One signal finding is that the overall 
size of a polymer of length N scales as N'^r, where Vp is a 
dimension dependent critical exponent. This is reflected 
in the behavior of various observables, such as the aver¬ 
age end-to-end distance, and the radius of gyration (Rg) 
mm- The scaling exponent is known to be sensitive to 
the excluded volume interaction applied to part or all of 
the monomers. Other geometrical constraints applied to 
the chain, such as cyclization (where the first monomer 
is connected to the last one), leading to a “ring” poly¬ 
mer, affect only the prefactor for these quantities. Ring 
polymers are commonly found in many biological sys¬ 
tems e.g., bacterial and mitochondrial genomes, as well 
as DNA plasmids used in many molecular biology exper¬ 
iments [3]. Recently, ring polymers were also studied in 
the context of a model for chromosome territories in the 
nucleus of eukaryotic cells [4]. 

The conformational fluctuations of some polymer mod¬ 
els can be analyzed using the theory of random walks 
(RW) [nil]. In particular, the fluctuations of the polymer 
size are of physical and biological importance. In the cur¬ 
rent paper we study the distribution of sample-averaged 
monomer-to-monomer distance of ring polymers, both 
for ideal, noninteracting polymers as well as for poly¬ 
mers with excluded volume. In particular, we exploit 
recent mathematical development on d dimensional con¬ 
strained Brownian motion (defined below) to find an 
exact expression for the distribution of sample averaged 
monomer to monomer distance for both ideal ring poly¬ 
mers and those with an additional excluded volume inter¬ 
action applied at a single point. This observable yields 
insight on the sample averaged fluctuations of polymer 
sizes. An important ingredient of this calculation is the 
identihcation of the appropriate boundary conditions for 
the underlying equation, a variant of the Feynman-Kac 


equation, which depends on both the dimensionality and 
the interaction, and in turn yields a selection rule for the 
solution. The resulting distributions are then compared 
numerically to those measured in simulations for a ring 
polymer with full excluded volume constraints. 

As we have noted, constrained random walks lie at the 
heart of our analysis. For rings, the primary constraint 
is that the path returns to the origin after N steps. 
Statistics of such constrained one dimensional Brownian 
paths have been the subject of much mathematical 
research [SU- These constrained paths have been given 
various names, depending on the additional constraints 
imposed. The basic case is that of a Brownian bridge 
where the return to the origin is the sole constraint. For 
a Brownian excursion, the path is also forbidden from 
reaching the origin in between. Majumdar and Comtet 
used Brownian excursions to determine statistical prop¬ 
erties of the fluctuating Edwards Wilkinson interface in 
an interval [muni. The focus of most previous work 
has been on the constrained one dimensional Brownian 
paths which describes inherently non-interacting systems 
(note that the problems of non-intersecting Brownian 
excursions [12] or vicious random walkers |13| in higher 
dimensions are exceptions). For the case of the fluctua¬ 
tions of ring polymers, we need to extend the theory of 
Brownian excursions and bridges to other dimensions. 
We address the influence of different kinds of interactions 
on the polymer structure, both analytically (for a single 
point interaction) and numerically (for a polymer with 
excluded volume interactions). These models yield rich 
physical behaviors and open new questions. 

Polymer Models. We consider three lattice models of 
ring polymers with N bonds, each of length b, in d di¬ 
mensions. The simplest polymer model is an “ideal ring” 
- a closed chain without excluded volume, where different 
monomers can occupy the same lattice site. While such 
a polymer does not exist in nature, its global behavior is 
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the same as that of a polymer at the Flory 0-temperature 
[I]. An ideal ring chain corresponds exactly to an unbi¬ 
ased RW in d dimensions which starts and ends at the 
origin, i.e. a d-dimensional bridge. In the second model, 
the “weakly interacting ring polymer”, the first and last 
monomers are tied to the origin and no other monomer is 
allowed to occupy this lattice site. This case is equivalent 
to that of d-dimensional excursions. The third model is 
a ring polymer with excluded volume interactions, also 
called a self-avoiding walk (SAW). Further details on the 
polymer models and simulation methods are provided in 
the supplementary material (SM). We first consider the 
ideal and weakly interacting ring models, for which we 
can provide an analytic solution. 

Bessel process. In the analogy between statistics of an 
ideal polymer and a RW, the position of the ith monomer, 
Ti, corresponds to the position of the random walker 
after i time steps and N is proportional to the total ob¬ 
servation time. The Bessel process [HI [15] describes the 
dynamics of the distance r = |r| from the origin of a 
Brownian particle in d dimensions. This process is de¬ 
scribed by the following Langevin equation: 

T- = ( 1 ) 

where ri(t) is Gaussian white noise satisfying (ry) = 0 and 
{'r]{t)ri{t')) = — One may map the polymer models 

to the Bessel process using (r(t)^) = dt = Nb"^ = (R^), 
where is the mean square end-to-end distance of an 
ideal linear polymer chain without constraints. In what 
follows we choose = d and t = N. 

Bessel Excursions and Reflected Bridges. The pro¬ 
cess r(t) with the additional constraint of starting and 
ending at the origin, is called a reflected (since r > 0) 
Bessel bridge. This process corresponds to an ideal (non¬ 
interacting) ring chain. Bessel excursions are paths still 
described by Eq. Q however with the additional con¬ 
straint that any path that reaches the origin (besides 
t = 0 and t = N) is excluded. The Bessel excursion cor¬ 
responds to what we have called the “weakly interacting” 
ring chain, where a multiple occupation of the origin is 
not allowed. The mapping of the polymer models to the 
Bessel process, allows us to extract statistical properties 
of the former with new tools developed in the stochastic 
community HISKII]. 

The Observable A. For a ring polymer, let be the 
position of the ith monomer in space where i = 0..N, and 
we place the origin at the position of the zeroth monomer, 
tq. For the weakly interacting chain this monomer is also 
the excluding one. We study a new measure. A, for the 
size of a ring polymer, defined by 

N 

A = ^|ri-ro|. (2) 

i=0 

In the RW language, A is the area under a random pro¬ 
cess, and hence is a random variable itself. Clearly I = 



A/(A) 


FIG. 1. (color online) Theoretical P+(A/{A)) Eq. (101 in di¬ 
mension d = 2,3,5 nicely matches simulations of the ideal ring 
(IR) and the weakly interacting (WI) models, the exception 
being the weakly interacting model in d = 2 due to critical 
slowing down (the simulations did not converge for N = 10®, 
see text). The theoretical curves (solid lines) for the two mod¬ 
els are identical for d > 2. The d = 1 WI theory is identical 
to the d — 3 theory, which the simulation confirms. The SAW 
simulations in d = 2, 3 are practically indistinguishable from 
each other and from the theoretical curve corresponding to a 
noninteracting ring in dimension d = 5. 


A/N is the sample averaged distance of the monomers 
from the origin. Specifically, let the area under the ran¬ 
dom Bessel curve r{t) be denoted by Ab = fo r(t)dt (the 
subscript B is for Bessel). More generally the mapping 
of the processes implies that in the limit of large IV the 
distribution of Ab/(Ab) is identical to the distribution 
of A/{A) (or I/{I)), with the corresponding constraints. 

Numerical results. In Fig. we plot the probability 
density function (PDF) of the scaled random variable 
A/(A). Both results of simulations and theory are pre¬ 
sented, however at this stage let us focus on the main 
features as revealed in the simulations. For dimensions 
d = 2,..., 5, there is a clear trend of narrowing of the 
PDF for increasing d. This trend is explained by examin¬ 
ing Eq. (1): As d increases, the noise term becomes neg¬ 
ligible compared with the force term, resulting in smaller 
fluctuations and narrower tails. Against this expected 
trend are the results in d = 1 for the weakly interacting 
model. As we will show analytically, the weakly interact¬ 
ing model in dimensions one and three surprisingly have 
the same distribution even though d = 1 has a vanishing 
deterministic force term in Eq. (1) while for d = 3 the 
force is clearly not zero. In addition, we observe that for 
d > 3 the shape of the distribution of weakly interact¬ 
ing and ideal ring chains coincide, indicating that weak 
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interactions are negligible (when N —>■ oo). As we shall 
see, this is also observed in the theory. Indeed the theory 
discussed below suggests that these two distributions are 
already identical for d = 2. However, since this is a crit¬ 
ical dimension, due to extremely slow convergence, we 
don’t see this behavior in the simulations. This asymp¬ 
totic convergence is logarithmic (see SM) and an £ expan¬ 
sion shows that it is reminiscent of critical slowing down. 
As for the SAW polymer, we see that the fluctuations 
are considerably reduced compared to the other models. 
This is due to the fact that the number of configurations 
of a SAW polymer is smaller than for the other models, 
hence fluctuations are smaller. A striking observation is 
that the two and three dimensional SAW results are iden¬ 
tical, both being equal to the simulations of the d = 5 
models. We now address these observations with theory. 

Functionals of Constrained Bessel Processes. Our goal 
is to find the PDF P{AB,t) of the functional Ab = 
/p r{t)dt of the Bessel process, constrained to start and 
end at the origin. We show that the difference be¬ 
tween the weakly interacting model (the Bessel excur¬ 
sion) and the ideal polymer (the reflected Bessel bridge) 
enters through the boundary condition in the Feynman- 
Kac type of equations describing these functionals. The 
choice of boundary condition turns out to be non-trivial 
and controls the solution. Other aspects of the solution 
follow the steps in [S]. 

It is useful to find first the Laplace transform of 
P{AB,t), i-e-, P{s,t) = P(AB,f)exp(-ABs)dAB to 
solve the equations, and invert P{s,t) back to P{AB,t)- 
Let Gt{r, AbIcq) be the joint PDF of the random pair 
(r. As) with initial condition Go(r, AbIcq) = 5{AB)5{r — 
ro) and G = Gt{r,s\ro) its Laplace pair. The modified 
Feynman-Kac equation reads m- 


1 

2 


dr'^ 


d l-d 
dr r 


G-srG=-G. 

at 


( 3 ) 


with G\t=o = Sir — ro) and tq a cutoff which is eventually 
taken to zero. For d = 1, the second term on the right 
hand side vanishes, and we get the celebrated Feynman- 
Kac equation corresponding to Brownian functionals m- 
The third linear term —srG stems from the choice of our 
observable, namely our functional is linear in r HH]. 
Since we are describing a ring polymer, the Bessel process 
must start and end on the origin, and so, following m, 
we need to calculate 


P(s,t) 


hm 

r=ro^O Gt(r,0|ro) 


( 4 ) 


The denominator gives the proper normalization condi¬ 
tion. 

The hrst step in the calculation is to perform a simi¬ 
larity transformation: 


G't(r,s|ro) = (^ —) Kt(r,slro). (5) 

Using Eq. (1^, Kt(r,slro) is the imaginary time propa¬ 
gator of a Schrodingier operator: 

ffKt(r, s|ro) -f ^.Kt{r, s|ro) = 5{r - ro)S{t) (6) 


with the effective Hamiltonian: 


y . 1 , ( d - 2)^-1 

2 dr"^ 


( 7 ) 


The effective Hamiltonian reveals a subtle symmetry, 
namely two systems in dimensions di and d 2 satisfying 
di + ^2 = 4 behave identically. Note that this symmetry 
is not affected by the choice of functional (or observable) 
since the latter only modifies the last term in H. This 
explains the identity of the d = 1 and d = 3 PDFs noted 
earlier. 

Boundary Conditions for Ideal and Weakly Interacting 
Models. The solution of Eq. ^ 

Kt{r, 5|ro) = ^ (8) 

k 


is constructed from the eigenfunctions (fk of H where 
Afe is the fcth eigenvalue and the normalization condition 
is /q°° <(>fe(?')dr = 1. The subtle point in the analysis is 
the assignment of the appropriate boundary condition 
corresponding to the underlying polymer models we con¬ 
sider. The eigenfunctions at small r exhibit one of two 
behaviors: 




dtr 


l+|2-d| 


or 






( 9 ) 


From the normalization condition, the </)“ solution can¬ 
not be valid for d > 4 and d < 0. For the critical dimen¬ 
sion d = 2 the two solutions are: ~ d^rs or 

dfr^ Inr. We now solve the problem for the two bound¬ 
ary conditions and then show how to choose the relevant 
one for the physical models under investigation. 

The distribution of A/{A). Following the Feynman- 
Kac formalism described above and performing the in¬ 
verse Laplace transform [3, we find two solutions for 
the PDF of the scaled variable y = A/(A) 
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( 10 ) 


The solution is independent of N and valid in the limit 
of —)■ oo. Here, |a| = \d — 2|/2, \v\ = 2|a|/3, and 
2 F 2 {-) refers to the generalized hypergeometric functions. 
The supplementary material provides a list of \k and dk 
values for d = 1,... 4. For d = 1 the solution agrees with 
the known results [ini El EOl El], where the + solution is 
the celebrated Airy distribution [iniE]- The average of 
A is 


(A)± = c±iV3/2, 


4v^r +1) 


( 11 ) 


The + solution was previously presented in a slightly 
different form in and here the question is how to 
choose the solution for the corresponding polymer mod¬ 
els. Clearly, for d = 2, (A)+ = (A)“, indicating that 
this is a critical dimension. Further (A)+ in 1 and 3 
dimensions are identical and so is (A)“ as the result of 
the symmetry around d = 2 in Eq. Q. The scaling 
(A) oc is expected since r scales with the square 

root of N as for Brownian motion, so the integral over 
the random processes r scales like 

We investigate the physical interpretation of the 
two possible boundary conditions. A mathematical 
classification of boundary conditions was provided in 
usmn] and here we find the physical situations where 
these conditions apply. We examine the behavior of the 
probability current associated with the mode: = 

(^’o) (f-))' + (1 - 

for r near the boundary ro —> 0 in dimension d. The 
analysis is summarized in Table [T| We see that in 
dimension two and higher, the current on the origin 
is either zero or positive. A positive current at the 
boundary means that probability is flowing into the 
system, which is an unphysical situation in our system. 
Hence we conclude that in dimension two and higher, the 
— solution is not relevant. This implies that statistics of 
excursion and reflected bridges (and equivalently, ideal 
and weakly interacting ring polymers) are identical for 
d > 2 and correspond to the -I- solution. 

In Fig. [^we compare the results of the simulations of 


1 

d 



Jt 


Jk 

d= 1 


^ r 

4+<o 

~ Const 
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II 

II 


Vr 

O 

II 

f’k ~ Vblogr 

Jk ^ ^ 

CO 

II 


^ r 
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II 

~ Const 



TABLE I. Behavior of the probability eigenfunctions and 
probability currents in the proximity of the origin (r = O'*") for 
the -f/— solutions in different dimensions. A current > 0 
on the origin is unphysical hence the critical dimension where 
local interactions are uniportant is 2. 


the ideal and weakly interacting polymer models with our 
theoretical results for P+(A/(A)), as given in Eq. (101. 
As noted above, for d > 3, we see that even for finite 
size chains the local interaction is not important, and 
that the theory and simulations perfectly match, while 
for d = 2 there are strong finite size effects in the weakly 
interacting case. 

Self-avoiding polymers. Extensive simulations of ring 
SAWs were performed on cubic lattices. As has already 
been pointed out, the global expansion of a polymer is 
characterized by the exponent Vp. Since A constitutes 
a measure of the overall size of a polymer, its behavior 
for large N should follow A ~ _/V'^p+i. For the ideal and 
weakly-interacting chains Vp = 1/2, and A For 

SAWs the exact value of the exponent depends on d, i.e., 
k'p = Vp{d). Vp = 1,0.75, and 0.5 are known exactly 
for d = 1,2, and 4 El], respectively, while for d = 3, 
zzp, based on renormalization group considerations and 
Monte Carlo simulations, is Vp ~ 0.588. These prediction 
were extensively tested numerically for the observable of 
interest A with a critical dimension of d = 4, character¬ 
ized by a very slow convergence of the weakly-interacting 
model (see SM). 

While the scaling behavior of the SAW model is differ¬ 
ent from that of the other two models (as reflected in t'p), 
as we have noted, the scaled PDFs, P{A/{A)) are never¬ 
theless similar. A striking observation is that the SAWs 
in d = 2 and 3 coincide to the precision of our measure¬ 
ments with the comparably narrow PDF of the d = 5 
non-interacting model, (see Fig. [^. That these distribu- 
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tion are narrower than the non-self-avoiding case can be 
qualitatively explained as follows: Since the interaction 
forbids many compact conformations, the fluctuations of 
the area become smaller. This is easily observed in the 
extreme case of a linear SAW in one dimension where 
only one conformation is allowed and the scaled PDF 
assumes the form of a 5-function. 

Discussion. The mapping of ring polymer models to 
the reflected Bessel bridge and excursion is very promis¬ 
ing since it implies that not only the observable A can be 
analytically computed, but also other measures of statis¬ 
tics of ring polymers. An example would be the max¬ 
imal distance from one of the monomers to any other 
monomer, since that would relate to extreme value statis¬ 
tics of a correlated process. The famed Airy distribution 
describes both the one dimensional polymer, as well as 
the three dimensional one, due to the symmetry we have 
found in the underlying Hamiltonian. The case of d = 2 is 
critical in the sense that interaction on the origin becomes 
negligible, though for finite size chains it is still impor¬ 
tant. Boundary conditions of the Feynman-Kac equation 
were related to physical models, which allowed as to se¬ 
lect the solutions relevant for physical models. The SAW 
polymer exhibits interesting behavior; the distribution 
of A/{A) is identical (up to numerical precision) in di¬ 
mension 2 and 3 and corresponds to the non-interacting 
models in dimension 5. Further work on this observation 
is required. 
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FIG. SI. (color online) {A)/N^^^ as a function of N for the critical dimension d = 2: weakly-interacting polymer simulations 
(blue circles) and the fitted curve of the function: {A)/N^^‘^ ~ 0.497 -I- 0.544/log(Af). The black dashed line is the theoretical 
value of c"*" = 0.4922 for N —^ oo. 


SUPPLEMENTARY MATERIAL 


Ring Polymer Simulations in d dimensions 


Ideal and Weakly Interacting Polymers 


In our simulations for the ideal and weakly-interacting ring polymers the chain is built of N consecutive bonds 
on a lattice. In the zth step, the bond displacement is = ±1 in each of the d directions, j = 1,..., d, hence 

the bond length is 6 = example, for d = 3, starting from the origin {i = 0) with first 


step of Ar = (-1-1,—1,-1-1) (yielding b = -s/S) we reach the lattice site ri = (-1-1,—1,-fl). For a second step of 
Ar = (—1, -1,-1) we end up at lattice site V 2 = (0, —2,0) for the i = 2 monomer. 

In order to maintain the closure condition of the chain, we choose an array of length of N with N/2 components of 
(-1-1) and N/2 of (—1) in each of the directions (d such arrays), and then shuffle them for each direction separately. For 
each i, the components of our d dimensional step are the ith values of these arrays. The sum of all of the displacements 
in each direction is then naturally zero so that the last monomer is always positioned at the origin. For the ideal 
chain model we built 10® such conformations while for the weakly interacting we threw away all the conformations 
that crossed the origin prior to the final monomer. For each of the conformations we calculated A = Eto 1^* ~ 
where Tq = 0, and plotted the distribution of this parameter. 

According to Eq. (11) in the paper, for —>■ cx) we have (A) = By this we can check the convergence 

of the simulations to the theory as a function of N. At the critical dimension, d = 2 this convergence becomes very 
slow. In Fig. SI we plot {A)/N^/'^ for different N values (in a logarithmic scale) in d = 2, where (A) = c+ = 0.4922 


is the theoretical value for N ^ oo. One can observe the very slow logarithmic convergence to this value. 


Self-Avoiding Ring Polymers 

Monte Carlo (MC) simulations have been applied to self-avoiding ring polymer models [1] on square, simple cubic, 
and d = 4 hyper-cubic lattices. The polymer consists of N monomers (and N bonds), where the first monomer 
is attached to the origin of the coordinate system on the lattice, and the Nth monomer is the nearest-neighbor to 
the origin. At step j of the MC process, monomer fc (1 < A: < A — 1) is selected at random (i.e., with probability 
1/{N — 1)) and the segment of m monomers following k (i.e., k + l,k + 2, ,k + m) become subject to change in the 
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A/{A) 


FIG. S2. (color online) Theoretical P{AI{A)) for reflected Bessel bridges and Bessel excursions for d = 4 (black solid line) along 
with simulations of d = 4 SAW for N = 200 (blue diamonds), N = 400 (green triangles), N = 800 (blue circles), N = 1600 
(purple squares). 


MC process; the rest of the chain (i.e., monomers 1 to fc and fc + to + 1 to N) is held fixed (notice that if k is at the 
end of the chain, — TO + l<fc<A^—1 ,to decreases correspondingly from to — 1 to 1). Thus, this current segment 
is temporarily removed and a scanning procedure is used to calculate all the possible segment configurations of to 
monomers satisfying the excluded volume interaction and the loop closure condition (i.e., the segment of to monomers 
should start at k and its last monomer, fc + to is a nearest neighbor to monomer fc + to + 1; notice that the initial 
segment configuration is generated as well). The segment configuration for step j is chosen at random out of the set 
of C configurations generated by the scanning procedure and the MC process continues. 

This process starts from a given ring configuration, whose transient influence is eliminated by a long initial simula¬ 
tion, which leads to typical equilibrium chain configurations. Then, every certain constant amount of MC steps the 
current ring configuration is stored in a file to create a final sample of n rings from which the averages and fluctuations 
of the physical properties of interest are calculated. The segment sizes used are to = 10 for the square lattice, to = 6 
for the simple cubic lattice, and to = 4 for d = 4. For each lattice several chain lengths, N are studied. 

We calculated the averages of Rg, 


N-l 




(SI) 


i=0 


where 


1 A^-l 

rc.m. = 

2=0 


(S2) 


ro|. For large N, (Rg) = ^ {R^) increases as ~ where Vp is a critical exponent and as 


and A = |r* - 

discussed in the text, i^a = + 1. However, we are mainly interested in the fluctuations of A, i.e., in the shape of 

the scaled distribution of A/(A) for different N. To check the reliability of the simulations we provide in Table SI the 
results obtained for Rg. The calculated results for d = 2 and d = 3 are equal within the error bars to those of Vpredicted 
while for d = 4 Vcaicuiated is too large due to a logarithmic correction to scaling, which would become insignificant 
only for much larger N. In fact, considering this correction in the analysis has led indeed to Vcaicuiated — 0.5. The 
same quality of results for va has been obtained for the observable A (where ~ 1.5 for d = 4). 

The d = 4 SAW case is a critical one, since the critical exponent, Vp for lower dimensions significantly differs from 
the i/p = 1/2 of the non-interacting models, and for d > 4 the interactions become unimportant for N —>■ oo. Hence, 
we expect P{A/{A)) of the d = 4 SAW to coincide with that of the non-interacting model. However, for finite N the 
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d 

N range 

^calculated 

^•predicted 

2 

200 - 1600 

0.7507 ± 0.002 0.75 (exact) 

3 

60 - 1002 

0.5890 ± 0.002 

0.588 

4 

200 - 2560 

0.525 ±0.015 

0.5 (exact) 


TABLE SI. 


Dimension Value 

fc = 1 

k^2 

CO 

II 

II 

k = 5 

II 

k = 7 

1,3 

K 

2.3381 

4.088 

5.2056 

6.7871 

7.944 

9.02265 

10.040 


dt 

1 

1 

1 

1 

1 

1 

1 

1 

K 

1.0188 

3.2482 

4.8201 

6.1633 

7.3722 

8.4885 

9.5345 


dk 

0.99088 

0.5550 0.4554 0.4241 

0.3837 

0.3663 

0.3566 

2 


1.738 

3.671 

5.170 

6.475 

7.658 

8.755 

9.787 


dt 

1.1391 

0.9195 

0.8386 0.7885 

0.7597 

0.7317 

0.7109 

4 

K 

2.873 

4.494 

5.868 

7.098 

8.231 

9.291 

10.294 


dt 

0.7585 

0.8807 0.9523 

1.0047 

1.0482 

1.0820 

1.1108 

5 


3.362 

4.885 

6.208 

7.406 

8.516 

9.558 

10.547 


dt 

0.5187 

0.6762 

0.7838 0.8646 

0.9405 

1.0012 

1.0531 


TABLE S2. The first 7 eigenvalues, A*, and the corresponding numeric coefficients dk required for plotting the theoretical PDFs. 
The eigenvalues of the + solutions for d = 1 and d = 3 are the negatives of the zeros of the Airy function: Ai(—A^) = 0. 
The eigenvalues A^ of the — solution are the negatives of the zeros of its derivative: Ai'(—A^) = 0. The d^ for d = 1 are 
tabulated in Ref. [^. 


interaction still has an effect on the distribution’s shape, and an even more pronounced one for ring polymers. For 
the values of N we used in our SAW simulations the curve had not yet converged as can be seen in Fig. (S2). A 
downward trend of the curves towards that of the non-interacting case (i.e. towards convergence) can nevertheless be 
seen. A similar problem is not found for SAW in d = 2,3 which are reported in the main text. 


Numeric values of Xk and dt 


The theoretical PDFs for different dimensions, presented in Eq. (10) in the paper, may be plotted using 
MATHEMATIC A®. In order to find the numerical coefficients Xk (eigenvalue) and dk (the normalization coeffi¬ 
cient of the eigenfunction) values of the fcth mode, we used the numerical method described in detail in [5] . In Table 
(S2) we present the values of the first few Xk and dk for different boundaries in different dimensions. We found that 
the first 7 eigenvalues are usually sufficient for the evaluation of P^{A/(A)). 
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